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Abstract 

We study different models of dark matter distribution for the halo of 
our galaxy. In particular, we consider Eddington and King-Michie 
models, which include anisotropy in the velocity space, and com- 
pute in a self-consistent way the amount of dark matter present 
in the halo. Assuming that the dark matter is in form of Massive 
Astrophysical Compact Halo Objects (MACHOs), we find for each 
model the expected number of microlensing events and their aver- 
age time duration for an experiment monitoring stars in the Large 
Magellanic Cloud (LMC). The main effect of including anisotropy 
is to reduce the microlensing rate by about 30% and to increase, 
but only slightly, the mean event duration, as compared to the stan- 
dard halo model. Consideration of different luminous models for the 
visible part of the galaxy also induce variations in the microlensing 
results by roughly the same amount as mentioned above. The main 
uncertainty, in order to be able to discriminate between different 
dark matter distributions and to estimate the fraction of it in form 
of MACHOs, is due to the poor knowledge of the rotation velocity 
at large galactocentric distances. 
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1. INTRODUCTION 



An important problem in astrophysics lies in the knowledge of the na- 
ture of the non-luminous matter present in galactic halos as well as of their 
shape and extension. The dark matter in galactic halos may be in the form of 
MACHOs (Massive Astrophysical Compact Halo Objects), in the mass range 
10"'^ < M/Mq <~ 10"^ (De Riijula et al. 1992) , which can be detected using 
the gravitational lensing effect (Paczynski 1986). The few microlensing events 
found so far by monitoring stars in the Large Magellanic Cloud (LMC) by the 
EROS (Aubourg et al. 1993) and the MACHO (Alcock et al. 1993) collabora- 
tions, while confirming the presence of MACHOs, do not yet allow to make a 
precise estimate of the fraction of the halo dark matter in the form of MACHOs 
nor to infer whether they are located in the halo or rather either in the LMC 
itself (Sahu 1994, Wu 1995, Gould et al. 1994) or in a thick disk of our galaxy. 
Assuming a standard spherical halo model, Alcock et al. (1995a) found that 
MACHOs contribute a fraction 0.19^^;^^ to the halo dark matter, whereas their 
average mass turns out to be ^ O.OSM© (Jetzer 1994). 

In this paper wc investigate different models of dark matter distribution for 
the halo of our galaxy. In particular we consider the so-called Eddington and 
King-Michie models, which include anisotropy effects in the velocity space. For 
these models we determine in a self-consistent way the amount of dark matter 
present in the halo. Moreover, assuming that the dark matter is in form of 
MACHOs we find for each model the expected number of microlensing events 
and their average time duration for an experiment monitoring stars in the LMC. 
These results are of relevance in order to determine the fraction of the halo dark 
mass in form of MACHOs. If the preliminary results of Alcock et al. (1995a) 
mentioned above will be confirmed by future improved observations, the problem 
arises of how to explain the nature of the remaining fraction of the halo dark 
matter. Besides the possibility of being composed of new exotic particles, it 
could still be baryonic and in form of molecular clouds (mainly H2). This latter 
scenario has been recently investigated in detail by us in several papers (De 
Paolis et al. 1995, 1995a) and also by Gerhard and Silk (1995). 

In section 2 we describe; the models for the dark matter distribution and in 
section 3 wc give the formiilas for computing the microlensing rate, the average 
time duration of an event and the average mass of the MACHOs using the mo- 
ment method. The numerical results and the conclusions are given in section 4 
and 5, respectively. 

2. MODELS FOR DARK MATTER DISTRIBUTION 

We assume that the galaxy contains two main components, namely, the 
visible component and the dark component. In the following we will assume 
that the dark matter consists of MACHOs, although the results on the total 
amount of the dynamical mass do not depend on this assumption and are. 
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therefore, of more general validity. In the last few years, the picture of the 
visible component of the Milky Way has evolved from the BSS model (Bahcall 
et al. 1983) which assumed stars to be distributed according to a central bulge, 
a spheroid and an exponential disk. In particular, the central concentration of 
stars is now described by a triaxial bulge model with the density law 

pc{x,y.z)^^e-^''\ with s^^ix^l-a^ + y^/hy + z^c^ (1) 

where the bulge mass is ^ 2 x 10"'^° Mq and the scale lengths are a = 1.49 
kpc, b = 0.58 kpc, c = 0.40 kpc (Dwek et al. 1994). The coordinates x and y 
span the galactic disk plane, whereas z is perpendicular to it. The remaining 
luminous matter (the spheroid and the disk of the BSS model) can be described 
with a double exponential disk (see Gilmore et al. 1989), so that the galactic 
disk has both a "thin" (Di) and a "thick" (D2) component. For the "thin" 
luminous disk we adopt the following density distribution 

p^(X,z) = ^e-N/^e-(^-«o)/^ (2) 

where the local projected mass density is Sq ~ 25 Mq pc^^, the scale parameters 
are H ~ 0.30 kpc and h ~ 3.5 kpc and i?o=8-5 kpc is the local galactocentric 
distance. Here X is the galactocentric distance in the plane. For the "thick" 
component we consider the same density law as in eq.(H) but with variable 
thicknesses in the range i? = 1 ± 0.5 kpc and local projected density Sq ~ 
50 ±25 Mq pc-2. 

The total local projected mass density within a distance of (0.3— 1.1) kpc of 
the galactic plane is measured to be in the range (40 — 85) Mq pc"'^. This ex- 
plains our chosen range of values for the "thick" disk, which in fact corresponds 
to a total luminous projected mass density of (50—100) Mq pc~^. In our models 
we also consider the effect of varying the bulge mass M^ = 2 ± 1 x 10^" Mq, the 
local galactocentric distance i?o — 8.5 ± 1 kpc and the scale length h ~ 3.5 ±0.5 
kpc, while we keep all the remaining parameters appearing in the previous equa- 
tions fixed. 

We treat MACHOs composing the dark matter with the formalism based on 
the equation of state assuming that they are spherically symmetric distributed 
and that: 

i) the velocities of MACHOs at any point of the galaxy are limited by the escape 
velocity v1{r) from the galaxy itself; 

ii) the distribution function can be anisotropic in the velocity space as a conse- 
quence of the initial conditions of the galaxy formation as well as of the MACHO 
formation processes [|. Then, for MACHOs of the same mass M, we adopt the 

^It is well known that for elliptical galaxies models in which the luminous matter is de- 
scribed by distribution functions with anisotropy in velocity space can account for their shapes. 
Moreover both isotropic and anisotropic models can explain optical and X-ray observations 
(De Paolis et al. 1995b). 
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King-Michie distribution function (King 1966, Michie 1963), which can be writ- 
ten as (Ingrosso et al. 1992) 

dn{r) = A(27ra2)-3/2 e^Wir)-WiO)] ^^-^^ I2a- _ g-WM) g-L^M/2L^ ^ (3) 

for V < Vc{r) and dn — otherwise. Here W{r) = w^(r)/2(T^ is the energy cutoff 
parameter, L = r x p the angular momentum and Lc the angular momentum 
cutoff. Moreover A is a normalization constant and a a parameter which in the 
limit of the classical statistics {W{r) oo and oo) represents the one- 

dimensional MACHO velocity dispersion. It is clear that the above distribution 
function gives lower values for the phase space density at large values of energy 
and angular momentum as compared to the ones obtained from the Boltzmann 
distribution. Furthermore it introduces an anisotropy in the velocity space, 
which increases with the radial coordinate r and leads to highly eccentric orbits 
for the MACHOs located in the outer regions of the galaxy. 

Visible and dark components are considered to be in hydrostatic equilibrium 
in the overall gravitational potential V solution of the Poisson equation W'^V = 
—AttG{ph + Pc + PDi + PD2)j which we solve assuming, as stated, spherical 
symmetry for the dark mass distribution. Clearly, this assumption as well as 
the purely radial anisotropy in the phase space for MACHOs are first order 
approximations, since the presence of the "thin" and "thick" stellar disks distort 
both the density and the velocity distribution of MACHOs close to the disk. 
However, we are interested mainly in the microlensing rate expected looking 
towards the LMC and we have verified that the errors in the number of events 
due to MACHOs close to the disks (i.e. in the region where the errors due 
to our approximations are higher) are negligible respect to those due to the 
uncertainties in the determination of the visible component. 

The distribution function given in eq.@ is the most natural way to build 
self-consistent galactic models, since the usual Boltzmann statistics requires an 
external cutoff (in the density) to avoid the mass divergence. Vice versa eq.(||) 
naturally implies a MACHO mass density equal to zero at the boundary R of 
the galaxy where W{R) = 0. Starting from eq.(||) and expressing Lc = Mvacr 
in terms of the anisotropy radius Va (Binney & Tremaine 1987), we obtain, after 
an integration over d^v, the radial component for the mass density 

Pff (r) = A (27ra2)-3/2 e[>^W->^(0)l (^) j^^'\e-i - e-^^^^^] F{\)d^ , (4) 

where A = {r/ra)^/£,, ^ = w^/2cr^, F{X) is the Dawson integral (Abramowitz & 
Stegun 1965) and W{r) = -V{r)/a^. 

The distribution function in eq. (^) can be approximated in the limit without 
energy cutoff {W — > oo) with a well treatable analytical model. The so-called 
Eddington model (Binney & Tremaine 1987) for which 
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In this case one can easily show that for values of r greater than the halo core 
radius a (the region in which the galactic dynamical mass is dominated by the 
halo dark matter) , the dark mass density pn (r) can be approximated by (Binney 
& Tremaine 1987) 

PH[r) = po 2 1 2 1 I / / ^2 ' (6) 

where po is the local dark mass density. 

Before going on with the presentation of the numerical results for the dark 
matter distributions, we discuss how we compute the microlensing rate and the 
average time duration. 

3. MICROLENSING RATES AND MASS MOMENTS 



When a MACHO of mass M is sufficiently close to the line of sight between 
us and a star in the LMC, the light from the source suffers a gravitational 
deflection and the original star brightness increases by 

u{v? + 4)1/2 V / 

Here u = cL/Re {d is the distance of the MACHO from the line of sight) and 
is the Einstein radius defined as: 

Re = - ^) ' (8) 

&■ 

with X — s/D, where D and s are the distance between the source (a star in 
the LMC), respectively the MACHO and the observer. 

Microlensing rates depend on the mass and velocity distribution of MA- 
CHOs. The mass density at a distance s from the observer is given by eq.(^ or 
(||). However, the adopted halo model does not determine the MACHO num- 
ber density as a function of mass. A simplifying assumption is to let the mass 
distribution be independent of the position in the galactic halo, i.e., we assume 
the following factorized form for the number density per unit mass dn/dM , 

with p = M/Mq, no{p) not depending on x and is subject to the normalization 
J dp^j^M = po- Nothing a priori is known about the distribution dno/dp. 

A different situation arises for the MACHOs velocity distribution, whose 
projection in the plane perpendicular to the line of sight f(x,VT) is obtained 
from eq.(^), by adopting cylindrical coordinates along the microlensing tube (so 
that d^v — VTdvTdvxdO). After integration over the longitudinal velocity Vx 
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and the direction 9 of the transverse velocity vt in the perpendicular plane we 
get 



(10) 

where 7 is the angle between the MACHO radial distance r and the distance s 
from the observer. The average value < > is obtained by a further integra- 
tion ^ ^ ^ 

< v^{i) >= [ f{5:,VT)v^^''+^UvT ■ (11) 



In the case of anisotropic models with W ^ 00 (the Eddington model in 
eq.(|^)) the previous equation can be cast into the following form: 

.p)„m /-oo 

< (i) >= T , , Z ■ 2 11/9 / Waii)] e-" 7y(™+i)d?? 



[l + (r/r,)2 sin27]i/2 

(12) 

where /q is a Bessel function and g{x) is given by 

{rlraf[l + {r/raf] sin^7 
[1 + {r/raY sm 7] 

In order to find the rate at which a single star is microlensed with amplifi- 
cation A > Amin (or u < Umax), wc considcr MACHOs with masses between ^ 
and II + dfj,, located at a distance from the observer in the range x and x + dx 
and with transverse velocity in the interval vt and vt + dvx ■ The collision time 
can be calculated using the well-known fact that the inverse of the collision time 
is the product of the MACHO number density, the microlensing cross-section 
and the average velocity. The rate dV at which a single star is microlensed in 
the interval dfidxdvT is given by (De Riijula et al. 1991, Griest 1991) 

dT{x,iJ,,VT) = 2DrEUmaxVTf{S:,VT)[fJ'X{l - x)]^^^ H {x)^^diidxdvT , (14) 
with 

..= (1^)''^ 3.2x10^ km. (15) 

One has to integrate the differential number of microlensing events, dN^v = 
Ni,tobsdT, over an appropriate range for fi, vt and i, in order to obtain the 
total number of microlensing events which can be compared with an experiment 
monitoring N^, stars during an observation time tohs and which is able to detect 
an amplification such that A > Amin- 
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The range of integration for is < wt < Vc{x) for the King-Michie models, 
while the upper limit goes to oo for the Eddington models. Moreover, x ranges 
between and 1 and fi varies in the interval where the mass distribution dno /dfi 
is not vanishing. 

However, each experiment has time thresholds Tmin and Tmax and only 
detects events with T^m ^ T < T^^ax^ and thus the integration range has to 
be such that T lies in this interval. The total number of micro-lcnsing events is 
then given by 

iv„ = /<i«„e(T-r„„oe(T_-r,, (le) 

where the integration is over the full range of dfidxdvT mentioned above. T 
is related in a complicated way to the integration variables, because of this no 
direct analytical integration in eq.(|l^) can be performed. 

In practice in order to evaluate eq.([l^) we define an efhciency function eo{fj,) 
which measures the fraction of the total number of microlensing events that 
meet the condition on T at a fixed MACHO mass and thus takes into account 
the 8 functions in eq.([l^). A more detailed analysis (De Riijula et al. 1991) 
shows that eo(/i) is, to a very good approximation, equal to unity for possible 
MACHO objects in the mass range of interest. We will therefore set eo{fj,) = 1 
in the following. 

From the experimental lensing data it is possible to extract information on 
the MACHO mass distribution (De Riijula et al. 1991). It proves most useful 
to average powers of the duration T, i.e., to construct "duration" moments. 
Consider an experiment that observes N^, stars during a time tobs, that has 
recorded a set of microlensing events, each one with a duration T. For each 
event we get a value for the dimensionless variable r, defined as follows: 

r^^T=^Ml-^)]V^ (17) 
te vt 

(here wq is the local rotational velocity) and we construct the n-moment of r 
from experimental data as: < r" >= Ylevents 

The theoretical expectations for these moments are calculated as follows: 

< r" >= j dN,, T^ = V Umax < > , (18) 

with m= (n + l)/2 and 

V =2N^tobsDrEVo , (19) 

l{m) = £ ^ di VTfii, vtMI - i)r'H{i) , (20) 

<//">= / dti ^fi^. (21) 
J dfj. 
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The main point of eq. (|18|) is that the experimental value of < r" > allows to 
determine the moments of the MACHO mass distribution, which can be used 
to reconstruct the mass function itself. The mean local density of MACHOs 
(number per cubic parsec) is related to < /i'^ >. The average local mass density 
of MACHOs is proportional to < /i^ > (solar masses per cubic parsec). The 
mean MACHO mass M is then given by 

where < > and < t^^ > are determined through the observed microlensing 
events. We may use < /i^ > to compute the fraction / of the dark matter 
density po that has been detected in the form of MACHOs. Indeed, we have 

. M© 1 1 s / 7.9 X 10-3 Mq pc-3 

/ = — < /i^ >w 126 < > pc^ — 

Po V Po 

Finally, the event duration T can be expressed in terms of the 7 functions defined 
in eq.(^o|) as follows: 

vo 7(1/2) ^ ^ 

In order to quantify the expected number of events it is convenient to take as 
an example a delta function distribution for the mass, i.e. dno/dfi = S{fj, — p,)/ fi. 
The rate of microlensing events with A > Amin (or u < Umax), is then 

^{Ajnm) =i^Umax ^ 2DrEUjnaxY7 7=/ dVT dx v'^ f {x , Vt)[x{1-x)]'^^'^ H {x) 

Mq Jo Jo 

(25) 

Assuming the standard halo model for the MACHO density distribution 
(which one obtains by taking the limit — > oo in eq.(|^)) Pff(r) = Poio,^ + 
Rl)/{r'^ + a^) with po = 7.9 x lO'^ M© pc'^, a = 5.5 kpc, using for the LMC 
distance Z? = 50 kpc and for the angle between the line of sight and the direction 
of the centre of the galaxy a = 82°, we obtain (De Riijula et al. 1991, Jetzer 
1991). 

f ^4^10-13 1^ -0 \( I \f PO \ I 



s V215 km/sy \^^D/kpc ) V^-S x 10-3 M© pc-^y ^M/Mq 

(26) 

For an experiment monitoring Ni, stars during an observation time tobs the total 
number of events with a amplification A > Amin is: Nev{Amin) = Ni,tobs^iAmin)- 
In Table 1 we show some values of Nev for the LMC, taking tobs = 1 year, 
= 10'^ stars and Amin = 1-34. 

4. RESULTS 



7 



We first discuss the Eddington model whose results coincide, in the isotropic 
limit, with those given in Table 1 for the galactic standard halo with flat rotation 
curve up to the LMC. We get a set of Eddington models by varying the values 
of the parameters a, po and Va- The amount of luminous matter, as discussed 
in section 2, is in the range between 4.8 x 10^" Af© ("minimum disk model") 
and 1.1 X 10-^^ Mq ("maximum disk model"). 

For each set of parameters we properly compute the resulting rotation curve 
considering the circular speed of the exponential disks (Binney & Tremaine 
1987) and the contribution due to the bulge and dark halo. A given model is 
physically acceptable only if it leads to a flat rotation curve up to the distance 
of the LMC (with a local rotational velocity vq — 215 ± 10 km s^^). In order 
to check this requirement we perform a X^j^mc test, which gives a measure of 
the flatness of the rotation curve in the region 5 kpc < r < 50 kpc. This way 
we find a range for the parameters for which the corresponding models satisfy 
the x^LMC < 1 test. 

In Table 2 we report the mean values of the allowed range for the parameters 
a, po and r^. The errors quoted in the Table are such that ~ 80% of the allowed 
models, which fulfill the X^lmc ^ test, have parameters whose values lie in 
the range defined by the errors around the mean quantity. Table 2 gives also 
the number of microlensing events N^v (as in Table 1 we consider an experiment 
monitoring 10^ stars during 1 year and for the MACHO we assume a mass of 
1O~^M0), the average event duration T, the ratio [7(0)77(1)], which is related 
to the mean MACHO mass Af, as well as the amounts of dark mass 
and Mff°* which are inside the distance to the LMC and 250 kpc, respectively. 
The latter value is of course arbitrary and should be regarded as an illustration, 
given also the fact that the true extent of the halo is unknown. 

As can be seen from Table 2, the main effect of the anisotropy in velocity 
space (second line) is to reduce the amount of M^f"* substantially. The dark 
mass Mj^^^'-^ also decreases with increasing anisotropy, however by less than 
20%. This is due to the fact that the anisotropy in the velocity space becomes 
particularly important starting from a distance r > (see eq.(^). 

The decrease of Af^*^*-^ for anisotropic models compared to isotropic ones 
implies a smaller number of microlensing events N^v , since the MACHOs orbits 
are now more eccentric. This reduction is also partially due to a shortening 
of the average transverse MACHO velocity < vt > with increasing anisotropy 
which, on the other hand, leads to an increase of the event duration. Finally, 
as one can see from the last column of Table 2, for the anisotropic models there 
is a decrease of the ratio [7(0)77(1)] leading (by eq.(|2^)) to a smaller average 
MACHO mass M. 

Till now we assumed a flat rotation curve in the range (5 — 50) kpc. Actually, 
the galactic rotation curve is well measured only in the range (5 — 20) kpc 
(Merrifield 1992). Thus we consider Eddington models for which we relax the 
condition of a flat rotation curve in the range (20 — 50) kpc. To select acceptable 
physical models we follow the strategy outlined, e.g., in Gates et al. (1995). We 
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require that: i) the local rotational velocity is vo = 215 ± 10 km/s; ii) the total 
variation in Vrot{r) in the range 5 kpc < r < 20 kpc is less than 14%; iii) the 
rotational velocity Vrot{LMC) at the LMC is in the range (150 — 307) km s~^. 

The results obtained for these Eddington models are given in Table 3, where 
the first row corresponds to isotropic models. A comparison between Tables 2 
and 3 shows that relaxing the requirement that the rotation curve has to be flat 
at distances larger than ~ 20 kpc does not change much the previous results. 
The main effect seems to be a larger range of variation for our results, although 
in Table 3 there is a trend towards a decrease of both Af^^^*^ and N^v The 
values in Table 3 also show the effect of the anisotropy which increases the 
event duration and decreases the [7(0)/7(l)] ratio. Notice that in both Table 2 
and 3 the "minimum disk model" leads to the higher value for the dark mass 
(columns 4 and 5), whereas the "maximum disk model" correlates with the 
smaller amount in the mentioned range around the mean quantity. 

Finally, we consider King-Michie models, whose results are reported in Table 
4. Here, again we adopt the same strategy as in Table 2, considering only models 
which lead to a flat rotation curve up to the LMC distance (performing again a 

X^LMC < 1 test). 

As before the numerical results depend on the assumed amount of luminous 
matter, so we consider, as an illustration, the two extreme models corresponding 
to the "minimum" (rows 1 and 4) and the "maximum" (rows 3 and 6) disk 
together with an intermediate model (rows 2 and 5). 

The values in Table 4 for the King-Michie models confirm the result that the 
anisotropy in phase space decreases N^v and M, while it increases T. However, 
the variation is not so marked as for the Eddington models in Table 2. Moreover, 
contrary to Table 2, there is no correlation between the values of Mj'^'-^ and 
Ngy . This can be understood by noting that increasing the amount of luminous 
matter (from rows 1 to 3 for the isotropic models, and from rows 4 to 6 for 
the anisotropic ones) the dark matter central density decreases while the core 
radius a increases. Therefore, the dark matter density is higher in the central 
region of the galaxy as well as in the outermost regions, whereas it decreases 
in the space in between the Solar system and the LMC. As a result of this the 
amoimt of microlensing events Nev will be smaller. 

The results in Table 4 show also an increase of the total dark matter with 
increasing luminous matter (dominated by the disk). The increase of dark 
matter with a corresponding increase of luminous matter can be understood by 
noting that the halo density as well as the velocity dispersion of the dark matter 
(randomly in orbit about the common center of mass) must be sufficiently high 
to avoid collapse into the disk. In the following paragraph we give a qualitative 
argument to show that this result is physically reasonable. 

Consider, first, a disk with a few dark particles distributed in a sphere. This 
may be thought of as an inhomogeneous sphere which must finally collapse to 
the disk over a sufficiently long period (due to gravitational instability). What 
will happen is that as the particle in the "sphere" orbits around the centre of 
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mass it will repeatedly cross the disk. While it is sufficiently close to the disk 
the particle will be attracted towards it and hence the orbit will be pulled into 
the plane. Now consider the sphere with only a very slight excess density in the 
equatorial disk. The inhomogeneity will still cause collapse into the disk but in 
a very much longer period. Clearly, as the density ratio of halo to excess matter 
tends to unity the sphere will be stabilized, whereas if it tends to zero it will be 
progressively de-stabilized. In other words, the time for collapse is a function 
of the density ratio. It tends to infinity as pu/ Pdisk ~* 1 and it tends to zero as 
Ph/ Pdisk 0. 

Now we turn our attention to the results in Table 4. The larger scatter of 
the values in the Table compared to the ones in Table 2 can be understood 
as being due to the variation of the amount of luminous matter as well as to 
the fact that the King-Michie models have an additional parameter (the central 
energy cutoff parameter Wo) with respect to the Eddington models (for which 
Wo — > oo). Due to the presence of this additional parameter it is not possible 
to see a clear correlation between the microlensing results and the degree of 
anisotropy in velocity space, since the latter effect can be smaller than the one 
induced by the variation of Wo. The effect of anisotropy would show up clearly 
once we compare models which all have a given value of Wq . 

The difficulty of getting precise microlensing results for the King-Michie 
models becomes even more evident if we adopt the same strategy to select phys- 
ical models, which lead to Table 3, that is by considering the conditions i) - iii). 
In Figure la we plot Nf,y as a function of the flatness degree [vrot{LMC) / vq] of 
the galactic rotation curve. In Figure lb, the same plot is given for the event 
duration T. In both figures the external and internal ellipses define the regions 
in which the models satisfy the condition of having a flat rotation curve only up 
to 20 kpc and to the LMC distance, respectively. It is clear that the knowledge 
of the LMC rotation velocity is crucial for determining the number of microlens- 
ing events, besides the ambiguity due to the uncertainty of the amount of the 
galactic luminous matter. 

5. CONCLUSIONS 

Several authors have studied the problem of determining the number of the 
expected microlensing events or, equivalently, the optical depth to microlensing 
by considering different models for the mass distribution, both luminous and 
dark in the galaxy (see, e.g. Alcock et al. 1995b, Kerins 1995, Kan-ya 1995, 
Evans & Jijina 1994, Evans 1994). Alcock et al. (1995b) considered "power-law 
models" for the halo and found that the microlensing rate can vary by as much 
as a factor 10 with respect to the value one gets for the standard halo. Kan-ya et 
al. (1995) analysed axisymmetric "power-law model" and studied the variation 
of the optical depth to microlensing and find that it can vary within a factor 2.5 
compared to the standard spherical halo model. Similar conclusions are reached 
by all other authors. 
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The main point in our present work is that we analyse another class of models 
(Eddington and King-Michie) which have anisotropy in the velocity space and 
we determine the corresponding parameters in a self-consistent way as described 
in section 2. We find that with the present knowledge of the various parameters 
the variation in the expected number of microlensing events is at least within 
30% from the value one gets for the standard halo model (flat rotation curve 
up to the LMC). This factor can even increase if one allows for less restrictive 
conditions. The typical time duration on the contrary seems to vary less, as 
well as the ratio [7(0)77(1)] which is related to the average MACHO mass. The 
determination of the latter with the moment method seems thus to be quite 
robust, in fact the ratio [7(0)/7(l)], as compared to the value for the standard 
halo model, varies by at most ^ ±30%. The main source of uncertainty, in 
order to be able to discriminate between different dark matter distributions and 
to estimate the fraction of it in form of MACHOs is due to the poor knowledge 
of the rotation velocity at large galactocentric distances up to the LMC. 
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MACHO mass in units of Mq 


Mean ii^; in km 


Mean microlensing time T 


N 


10-^ 


3 X 10** 


23.5 days 


5.6 


10-2 


iqs 


7.4 days 


17.8 


10-^ 


10^ 


17.8 hours 


178 


10-6 


10^ 


1.8 hours 


1784 



a 


Po/10-^ 






Mf/°' 




T 


[7(0)/7(l)] 


(kpc) 


[Mq pC-3) 


(kpc) 


(10" Mq) 


(10" Mq) 




(days) 




6.3 ± 1.4 


7.6 ± 1.5 


oo 


4.9 ± 1.0 


26 ±5 


5.4 ±0.9 


22.4 ±0.7 


7.3 ±0.2 


6.5 ± 1.4 


8.3 ± 1.3 


47± 12 


3.7 ±0.8 


7±2 


4.1 ±0.6 


27.0 ± 1.1 


5.2 ±0.3 



a 

(kpc) 


{Mq pC-3) 


ra 
(kpc) 


(10" Mq) 


Mil°* 

(10" Mq) 




T 

(days) 


[7(0)/7(l)] 


6.2 ± 1.4 


7.0 ±3.8 


oo 


4.5 ± 1.5 


24 ±9 


4.7 ± 1.5 


22.5 ±0.9 


7.7 ±0.2 


6.4 ± 1.4 


7.6 ± 4.0 


40± 15 


2.9 ± 1.5 


5±2 


3.5 ± 1.3 


27.7 ± 1.5 


5.7 ±0.3 
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a 


Po/10^^ 


V 




Mf/'"* 


Nev 


T 


[7(0)/7(l)] 


(kpc) 


(Mq pC-3) 


(kpc) 


(10" Mq) 


(10" Mq) 




(days) 




4.9 ±0.9 


9.1 ±0.9 


oo 


6.0 ±0.9 


24 ±9 


6.4 ±0.9 


22.6 ± 1.0 


8.2 ±0.8 


6.6 ± 1.9 


6.4 ±0.7 


oo 


6.3 ±0.9 


32 ±9 


5.9 ±0.9 


21.7± 1.0 


8.5 ± 1.0 


13 ±4 


4.0 ±0.5 


oo 


6.9 ± 1.0 


49 ±9 


5.8 ±0.9 


19.2 ± 1.0 


10.2 ± 1.1 


5.7± 1.5 


9.7± 1.1 


40± 15 


4.9 ±0.7 


11 ±3 


5.7 ±0.8 


25.1 ±0.9 


6.9 ±0.6 


8.9 ±3.5 


7.0 ± 1.0 


42± 16 


5.1 ±0.8 


14 ±3 


5.2 ±0.8 


24.1 ± 1.3 


7.3± 1.1 


17±3 


4.5 ±0.6 


45± 13 


5.8 ±0.7 


21 ±3 


5.0 ±0.9 


21.6± 1.0 


8.6 ±0.8 
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FIGURE CAPTIONS 

Figure la: The number of expected microlcnsing events N^v is given as a func- 
tion of the ratio [vrot{LMC)/vo] for the King-Michie models. We consider an 
experiment monitoring 10^ stars in the LMC during an observation time of 1 
year and we assume, as an indication, a MACHO mass of 1G~^ Ad^. Models 
with a flat rotation curve up to the LMC lie inside the inner ellipse (see also Ta- 
ble 4), whereas the ones with flat rotation curve only up to 20 kpc are bounded 
by the outer ellipse. 

Figure lb: The same as in Figure la for the average event duration T (in days). 
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TABLE CAPTIONS 



Table 1: Galactic standard halo model. 

Table 2: Mean values of parameters for Eddington models with flat rotation 
curves up to the LMC (the different parameters are described in the text). 

Table 3: Mean values of parameters for Eddington models with flat rotation 
curves up to 20 kpc. 

Table 4: Mean values of parameters for King-Michie models with flat rotation 
curve up to the LMC. 
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This figure "figl-l.png" is available in "png" format from: 
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